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Abstract The properties of a Lanthanum bromide (LaBr3) detector and its response functions were investigated via 


experiments and simulations in this paper. The LaBr3 detector had good relative energy resolution and higher 


efficiency than a high-purity germanium detector. Monte Carlo and other numerical methods were used to calculate 


the efficiencies of a LaBr3 detector with a square collimation window. A model of the numerical method was 


established based on a pure geometric model that was consistent with the experimental situation. The results showed 


that the detector response functions calculated by these methods were in great agreement with experimental results. 
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1 Introduction 


The radionuclide activity in radioactive waste drums 
must be characterized before handling and disposal of 
them. Segmented gamma scanning (SGS) and 
tomographic gamma scanning (TGS) are two widely 
used non-destructive radioactive characterization 
methods. The accuracy of SGS method is lower than 
that of TGS method when the waste drum is filled with 
heterogeneous radioactive materials. However, SGS 
method consumes less time in measuring a radioactive 
waste drum compared with TGS"'!, At present, the 
long measurement time is the biggest constraint on 
applying TGS in nuclear power plants. Two 
approaches have been used to reduce the measurement 
time in TGS—one approach is to establish a new 
reconstruction algorithm and the other is to use a 
detector with higher efficiency. Liu used dynamic 
grids in TGS to reconstruct source distribution”. The 
measurement time was reduced by half, while remains 


the measurement accuracy. To detect prompt gamma 


Lanthanum bromide 
scintillator (LaBr3) was applied to PGNAA"!, In the 


present study, a lanthanum bromide (LaBr3) detector 


rays from neutro capture, 


was applied to TGS. Its efficiency and energy 
resolution are described in Section 2. Section 3 
discusses the efficiencies of the detector determined by 
Monte Carlo and deterministic calculation methods. 


The conclusions are given in Section 4. 


2 Properties of LaBr; crystalline and its 
detector 


2.1 Properties of LaBr; crystalline 


Lanthanum bromide scintillator (LaBr3) with energy 
resolution of approximately 3% at 662 keV provides a 
iodide 
scintillators in resolution". Unlike HPGe detectors 


substantial improvement over sodium 
which require low-temperature working conditions, 
LaBr; detectors can be operated at room temperature. 
However, LaBr3; has several disadvantages such as 
internal radioactivity and low-energy response. The 


internal radioactivity is due to naturally existing 
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88ra and decay products of Ach. 


radioisotopes 
Here, the following properties of a LaBr3 detector 
were discussed, including Energy linearity, Relative 
energy resolution, detection efficiency as a function of 


y-ray energy. 


2.2 Properties of LaBr3; Detector 


The LaBr3 crystal was manufactured by Saint-Gobain. 
The following a set of sources was used in the 
experiment: 'S3Ba, ‘Bu, Cs, Co, Ra and “Am. 
These sources were located at 25 cm away from the 
detector center, and the measurement time was 3600 s. 
These sources supplied homogenous and wide range 
energy of y-rays, which were used to analyze the 
ability of LaBr; detectors and distinguish between the 
different radionuclides. The experiments primarily 
aimed to analyze the absolute efficiency, energy 
linearity, and peak width at half height of the LaBr; 
detector. Fig.1 shows a typical spectrum with '°’Cs 
(661.7 keV) and Co (1173.21 keV and 1332.47 keV) 
measured by the LaBr3 detector. Both radioactive 
sources could be identified, however, the spectrum 
was affected by the self-activity of the detector (1430 
keV from '**La and 1465 keV from °K) and the 789 


138 


keV peak from ` Ce was determined by f-decay. 


Fig.2 shows a good linearity of the energy from 41 Am 
(59.54 keV) to '”Eu (1408 keV). This energy range 
covered the most important radioactive sources in the 
radioactive waste drums. Fig.3 shows the relative 
energy resolution from 24 Am (59.54 keV) to Co 


(1332 keV), which was fitted by the function: 


FWHM _494 0354 |1 +0.00161 (1) 
E, E, E, 


The relative energy resolution is improved 
from 13% at 59.54 keV to 2% at 1332 keV, which is a 
satisfactory value for scintillation detectors. Fig.4 


shows the full absorption efficiency of the LaBr3 and 
HPGe detectors from 121 to 1332 keV. Co supplied 
high y energy, whereas '3Ba and '°’Cs supplied low 
and middle y energies, respectively. The measurement 
time was 3 600 s for both detectors, and the count for 
each peak reached 10 000. The statistical error for 
each peak was less than 1%. The full absorption 
efficiency was fitted by Eq.(2). 


ae ETE 
E E (2) 
+3.258x107 E -1.607 °F? 


The ratio of the efficiency of the LaBr3 detector 
to that of the HPGe detector ranged from 1.3 to 2 as 
energy increased from 121 to 700 keV. The calculated 
ratio was 2 when energy exceeded 700 keV. 
Consequently, the measurement time of the LaBr3; 
detector was the half of that of the HPGe detector. 
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Fig.1 Typical spectrum recorded with a LaBr; detector using 
'7Cs and Co 
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Fig.2 Channel vs. energy data from 59 to 1408 keV for the 
LaBr; detector. The solid line denotes a linear fit to the data. 
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Fig.3 Relative energy resolutions for the LaBr3 detector from 
59 to 1332 keV. The solid line represents the fit to the data. 
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3 Experimental and simulation for the 
LaBr; detector response 


3.1 Experimental description 


In the experiment, the detector was collimated by a 
square window. The source was distributed on a 9-cm 
radius circle (r,) in front of the detector. The distance 
(do) between the circular center and the detector was 
52 cm. Fig.5 schematically illustrates the experimental 
setup of the collimated detector system. The crystal is 
1.96 cm in radius (ra) and 3.91 cm in length. The 
collimator is 12 cm in length (koi), and its square 
window (dco!) is 4.8 cm in width. The drum was 
rotated in steps. The source located at 9 and 20 cm 
radial positions to a location shifted 3.5 cm (/) to the 
drum center. A photon count rate was measured and 
the drum was rotated for 15° in each step. 
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Fig.4 Efficiencies for LaBr; and HPGe detectors from 59 to 
1332 keV. The solid and dashed lines denote the fit to the data. 


Fig.5 Geometric model of the collimated detector system of 
the experiment. |: the offset between the center of the source 
circle and the axis of the detector; A: the projection of the 
center of the source circle on the axis of the detector; B: the 
projection of the point source C on the axis of the detector; dx: 


the distance between the point source C and detector surface; dy: 


the distance between the point source C and the axis of the 
detector; d,: the distance between the point source C and 
detector surface center; dy: the distance between point A and the 
detector surface; and ø: the angle of the point source. 


3.2 MCNP simulation of the LaBr3 detector 


response 


This section presents the results of the simulations of 
the LaBr3 detector response using the Monte Carlo 
N-particle (MCNP) method. The efficiencies of the 
sources in different positions are basic parameters in 
gamma scanning techniques. It was difficult to 
determine the efficiencies via experiments because the 
energy from experimental sources cannot cover all 
energies in the waste drum. By contrast, the energies 
and positions of the point sources as well as the void 
sources can be easily established in the simulation 
Thus, 


simulations to a certain extent. 


model. experiments can be replaced by 

In Fig.6, the simulated values are compared 
with the experimental results. A good agreement 
between the 


efficiencies is obtained. The average error is less than 


experimental and the simulated 


3%, indicating that the simulation model can fit the 
experimental results very well. 
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Fig.6 Comparison of the experimental results (circle points) 
with the simulation efficiencies by MCNP (square points) for 
137Cs point source (E, = 0.661 MeV). 

3.3 Simulation of the LaBr; detector response by 


the numerical method 


3.3.1 
Although MCNP simulation can fit the experimental 


Numerical method 


results, it needs a long time to calculate the detector 
method 
significantly decreases the calculation time. For 


efficiency. By contrast, the numerical 
example, it takes about 9 minutes to calculate the 
efficiency of one point by MCNP while only 5 minutes 
to calculate the efficiency of 90 points by numerical 
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method using the same computer. Thus, the numerical 
method may replace the MCNP. A calculation model 
for a point source is shown in Fig.5. The detector 
efficiency can be calculated based on the solid angle 
subtended by the collimated detector to the point 
source. Fig.5 presents geometry in accordance with 
Krings and Mauerhofer (2011)'°! as well as Nan et al. 
(2012)!7! The absolute detector efficiency e(d,) to the 
point source can be calculated by 


dy 


ad) =atd) | § Sin. Cop (3) 


x det 


where e(d,„) is the absolute detector efficiency to a 
point source located at the distance d, to the detector; 
é(do) is the absolute detector efficiency of a point 
source located at the drum center; San is the part of the 
active detector surface illuminated by the photon beam 
due to the collimation; Saet is the area of the active 
detector surface; C is correction factor of edge 
penetration effect; m is the power of (do/d,). 

The efficiencies were calculated based on the 
efficiency of the source located at Point A. First, the 
efficiency of point A was transformed to Point B 
located at d, (Fig.5). This transformation mainly 
depended on the distance between the point source and 
the detector surface. Then, the efficiency of Point B 
was transformed to the target location (Point C). This 
transformation mainly depended on the active detector 
surface illuminated by the photon beam San. The factor 
Cor corrected the absolute efficiency considering the 
edge penetration effect. Unscattered part of the photon 
beam that passed through the collimator and reached 
the active detector volume were disregarded because 
they had negligible effects on the transformation of the 
efficiencies. The source was divided into large point 
sources when the source was a volume source. The 
efficiency of this volume source e represents the 
average efficiencies of the point sources. The 
efficiency e is calculated by 


1 i) Si 
e=—Ya(d,)-|“2| 2. Cor 4 
por o) (2) 5, or (4) 


where k is the number of point sources in the volume 


source. 


3.3.2 Illuminated active detector surface 
In this section, rg is the radius of the active detector 
surface. Fig.7 presents six cases of the schematic 
calculation of San referring to Nan et al. (2012)! 
Before calculating Sg, in these cases, the arc segment 
areas (S) must be calculated. The area of an arc 
segment according to the notation used in Fig.5 is 
given by 
2 

5-42 Se0s21,, (5) 
where g is the length of the chord given by the 
intersection of the active detector surface and one side 
of the collimated photon beam. The central angle of 
the sector is expressed as 

5 = 2arcsin £- (6) 

2r, 

The side length of the photon beam at the plane 
defined by the active detector surface is given by 


l 1 
a. =—“—_ xa 7 
pro d, _ = col ( ) 


where d, is the projection of d, on the X-axis. 
According to Eqs. 5-7 and Fig.7, Sa, for these cases 
can be calculated by 


0 case 1 

nr? case 2 

Si tila, =b lla, -b| case 3 (8) 
Sun = 2 i i i 

nr? —(S, + S,) case 4 


Fig.7 Active detector surface illuminated by the photon beam 
using a square collimation. S;: arc segment areas in different 
cases. 


3.3.3 Edge penetration correction 
Fig.8 schematically presents different path lengths of 
photons in a detector crystal from a point source. The 
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correction factor was calculated as the average 
absorption probability of the photon beam with respect 
to the average absorption probability of a photon beam 
detector. The 
calculation of the average absorption probability relied 
distributed 
homogenously at the active detector surface. The 


from a point source facing the 


on a mesh of no virtual points 
penetration lengths /;r, of photons from a point source 
were calculated for each photon beam that entered the 
active detector surface at one of the n points within the 
illuminated part of the active detector surface. The 


correction factor can be expressed as 


-1 
Cor = i Ase Jf -18e (9) 
no Ny ial 
n: the number of virtual points distributed at the active 
detector surface in different conditions; u: the linear 
attenuation coefficient of Ge; p: the density of the Ge; 
no: the maximal number of virtual points distributed at 
the active detector surface. 


Point source 


Fig.8 Scheme of the different penetration lengths of photons 
from a point source in a collimated detector crystal. 


3.3.4 


The detector efficiency ¢ is determined by the detector 


Parameter m 


intrinsic efficiency ége and the detector spatial 
efficiency Espa 

E= Eget X Espa (10) 

The detector 

proportional to the solid angle of the cone formed by 


spatial efficiency Espa is 


the point source and the detector surface. 
1 0 
Eg, =—| 1-cos— 11 
w=5 [ 2) (11) 


where @ represents the solid angle of the cone formed 
by the point source and the detector surface. 


According to Taylor series expansion, 


coal) mala) la) 
4\2) 8x4!l2 2 


when d, >> rg, 0 is expressed by 


Ox 2ra (13) 
2 
Egt 1 
ein. 14 
4 d (4) 
Ea (2) (15) 
Eg d, 


The distance between the source and the detector 
surface is always larger than the radius of the detector. 
According to Eqs. 3, m was 


l ) 

EB 
l ) 
d. 


Thus, the value of m equalled —2 


m= 


3.3.5 


The count rate distribution of '*’Cs point sources was 


Reliability analysis of numerical method 


simulated under the same experimental condition. The 
experimental and calculated results about the count 
rate are shown in Fig.9. The calculated results were in 
a good agreement with the experimental results for all 
simulated source positions. The average error was less 
than 3%. This verifies the rationale of the numerical 
method in simulation of the LaBr; detector response. 
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Fig.9 Comparison between the experimental results (symbols) 
to calculation efficiencies according to the model (line) for 
137Cs point source (E, = 0.661 MeV). 
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4 Conclusion 2. 
This study demonstrated that the LaBr3 detector can be 3. 
used in SGS and TGS techniques. The error between 
the experimental results and the MCNP simulation j 
results in the detector efficiencies was less than 3%. 
The numerical method used to calculate the response 5 
of the LaBr; detector geometry of the gamma scanning 
systems agreed perfectly with the experimental data. 
Therefore, the numerical method of the LaBr; detector 6 
response can be used in SGS and TGS techniques. a 
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